Magnetic field-induced gapless state in multiband superconductors. 
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We investigate theoretically the properties of s-wave multiband superconductors in the weak cou- 
pling (BCS) limit in the presence of pair-breaking effects of magnetic field, ft is shown that a 
qualitatively new gapless superconducting state must appear in quasi-2D superconductors in mag- 
netic fields parallel to the plane, corresponding to a Sarma state induced on one of the Fermi 
surfaces. The emergence of the new state in s-wave multiband superconductors in the absence of 
anisotropy or spin-orbit interaction is usually accompanied by a zero-temperature first-order meta- 
magnetic phase transition. For anisotropic or non s-wave multiband superconductors the order of 
the zero-temperature metamagnetic transition depends on model parameters, and it may take the 
form of a smooth crossover. The details of the temperature-magnetic field phase diagram for multi- 
band superconductors are investigated analytically at zero temperature and numerically at a finite 
temperature. It is shown the zero-temperatures first-order phase transition gives rise to a critical 
region on the B — T phase diagram. We suggest possible experiments to detect the new gapless 
state. 

PACS numbers: 74.20.-z, 74.25. Ha, 74.70.Tx, 74.70.Kn 



I. INTRODUCTION 

Superconductivity in multiband metals was first inves- 
tigated theoretically shortly after the BCS theory[l|, 0, 
S 0, S] ■ A resurgence of interest in multiband supercon- 
ductivity has been mostly related to experimental discov- 
ery of MgB2 , and observation of the two s-wave gaps 
by various techniques 0, i, i, 0, Ei E [H, [ll ■ Re- 
cent discovery of iron oxypnictides[l5J, a new family of 
quasi-2D high temperature superconductors, has also at- 
tracted enormous theoretical interest to the problem of 
multiband superconductivity. First principles numerical 
band structure calculations [la. IitI. [l8l| show that several 
bands cross the Fermi level in these materials. While the 
pairing mechanism in iron oxypnictides is not yet clear, 
several bands are involved in the determination of both 
normal and superconducting properties |l8l ll9| . The ex- 
istence of multiband energy spectrum '20l| and associated 
several energy g:apS}21i. _22] was also verified experimen- 
tally in an older, quasi-2D s-wave superconductor NbSe2. 
BCS investigations of multiband superconductivity were 
restarted[23, ^2J|j but were mostly centered around the 
new physics that arises due to the presence of two sep- 
arate gaps. The multiband energy spectrum is also 
present in most unconventional heavy fermion supercon- 
ductors due to extreme complexity of the band structure 
of these materials, as indicated by de Haas van Alfven 
measurements. For example , gapless superconductivity 
of Abrikosov-Gor'kov type|25J was recently observed in 
thermal conductivity data for La-doped CeCoIn5[26], as 
indicated by unusual Wiedeman- Franz 1/x behavior of 
thermal conductivity at small concentration of La, x. 
The unusual behavior was attributed to the multiband 
structure J26ll27l| of this material, a d-wave superconduc- 
tor. 

Theoretical study of multiband superconductivity has 
not only been motivated by the above compounds. A 
somewhat modified multiband model is applicable to 



other new materials [28|, where a single Fermi surface gets 
spin-split into two pieces with different Fermi momenta 
due to interactions, such as superconductors without cen- 
ter of inversion (CI) [23] (for example, CePt.^Si,3Q,|), or 
ferromagnetic superconductors (UGe2 [3l|, [33] , ZrZn2[33], 
or URhGe[3J]). A particular case of CLsymmetry break- 
ing is two-dimensional surface superconductivity }3Es. i36l] . 
where two Fermi surfaces arise as a result of spin-orbit 
interaction of Rashba formp37l.[38|. 

The Fermi surface will also get spin-split in exchange 
field, an external magnetic field without orbital effects. 
Theoretically, paramagnetic pair-breaking by exchange 
field corresponds to the old problem of unbalanced pair- 
ing, first studied in the early eO-sfH EO, EE El, El, El , 
when the B — T phase diagram for three-dimensional ma- 
terials in exchange field has first been obtained. Orbital 
effects are almost always present in three dimensions, 
and even in very anisotropic quasi two-dimensional (2D) 
materials, unless Hc2 in magnetic fields parallel to the 
2D plane close to Clogston paramagnetic limit [39l [40|. 
This is, perhaps, the reason why the consequences of 
the theory have not studied in detail experimentally. A 
study of thin films in magnetic fields parallel to the sur- 
face remains, perhaps, the most promising experimen- 
tal setup for the observation of unbalanced pairing in 
superconductors |45l [46l [47J . A possible observation of 
the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) [41 , 42] has 
also been reported in some quasi 2D superconductors, 
such as those based on charge-transfer organic salts of 
BEDT-TTF or ET - ion[48] and the 1-1-5 family heavy 
fermion materials [49|, [50|, [5l|, [5^, [53| . 

While most theoretical studies of superconductivity in 
multiband compounds concentrated on the mechanism of 
superconductivity and the effects of several energy gaps 
on superconducting properties, some recent studies |54l] 
found that a new class of superfluids could potentially 
arise in these materials, one that features coexistence of 
fully gapped and gapless states. For the most part, this 



class of states was proposed in Bose-Einstein condensa- 
tion for different non-identical ferniions condensed by an 
optical trap [531, or in high-energy physics [50|, where the 
"breached" superfluid state appears as a result of pair- 
ing between different quarks [5 lliSS] in the asymptotically 
free limit. 

Gapless solutions of this type were previously known, 
but in most cases were found to be energetically unsta- 
ble. The first example corresponds to the second unsta- 
ble solution in the unbalanced pairing problem, and is 
commonly referred to as Sarma state[43, |4J|- The sec- 
ond solution of the gap equation for a superconductor 
placed in an exchange field was first obtained by Green's 
functions method by Baltenspergerj59] and Gor'kov and 
Rusinov[60|, where it was found to be energetically un- 
stable. The exact nature of this unstable state, and 
the existence of unpaired electron- and hole- Fermi sur- 
faces was later clarified within ECS theory by Sarma 
and Takada[43l |44|). The second example is pairing in 
a doped (unbalanced) excitonic gas, leading to an ex- 
citonic condensate |6ll. |63 |. which is actually the same 
two-band unbalanced pairing problem as one considered 
by Liu and Wilczek[5a|. The solution previously found 
for this problem takes the form of a magnetically (ferro- 
magnetically or antifcrromagnetically) polarized gapless 
state|63, [6J], and is different from "interior gap super- 
fluid" of Liu and Wilczek. Similar to Sarma state, this 
solution is usually unstable with respect to formation of 
LOFF state or domains |65l. l66l. l67| . although the stability 
of this magnetic solution has not been fully investigated 
as a function of masses involved. The third example is 
the Sarma gapless state, which appears in the problem 
of s-wave pairing in ferro mag netic superconductors ^68J . 
While it has been claimed[68[ that Sarma state is stabi- 
lized by the presence of ferromagnetic order, this claim 
has been debated later[69;, iTQ]. According to recent work 
of Liu and Wilczek, a large difference in effective masses 
on two Fermi surfaces tends to stabilize the "breached" 
superfluid state already within the unbalanced pairing 
problem [53]. 

The main motivation for this paper is a detailed theo- 
retical study of possible new gapless states in multiband 
superconductors in the presence of exchange magnetic 
fleld|7l|, or the multiband unbalanced pairing problem. 
While certain similarity does exist, the multiband prob- 
lem is really different [T^] from the unbalanced pairing 
problem considered in the above cases; pairing between 
two different species of fermions (two different bands) 
is usually not a relevant mechanism in multiband su- 
perconductors, since the energy difference for the two 
bands, Ae ~ leV is much greater than Tc ^ IK. Nev- 
ertheless, for a superconductor placed in an external 
"exchange" field the unbalanced pairing problem is re- 
covered. Surprisingly, as it was shown in Ref.[7l|], the 
multiband structure often leads to a a stabilization of 
unusual Sarma state on the second band in exchange 
fields ^j-bB > A2, where A2 is the energy band on the 
second Fermi surface. Thus, in quasi-2D superconduc- 



tors where several bands cross the Fermi surface, such 
as CeCoIus of 2H-Nb5'e2, the pecuharities of the B — T 
phase diagram may not limited to high magnetic fields, 
where the LOFF-related phenomena are observed. New 
singularities and gapless states associated with the gaps 
on secondary Fermi surfaces must arise in low magnetic 
fields as wellTT]. They correspond to the appearance of 
Sarma statei4^^] on the Fermi surface with smaller gap, 
one that becomes energetically stable due to the pres- 
ence of the superconducting gap on the other Fermi sur- 
face. This state is characterized j43j by the presence of 
unpaired spin-polarized electrons near the Fermi surface 
of the second band, two open electron- and hole- Fermi 
surfaces, a paramagnetic magnetic moment, and a first- 
order phase transition that always accompanies the ap- 
pearance of this unusual state in s-wave multiband su- 
perconductors in low magnetic fields [7l|. 

The paper is organized as follows. In Section II we in- 
troduce the multiband model and generalize its known 
solution in the s-wave case to include effects of gap 
anisotropy and non s-wave pairing symmetry. We demon- 
strate that the effective coupling constants can be elim- 
inated in favor of the measurable parameters for multi- 
band superconductors, such as the superconducting tran- 
sition temperature, Tc, and the ratio of the gap ampli- 
tudes and the densities of states on different Fermi sur- 
faces. In particular, thermodynamics of multiband su- 
perconductors is additive; the thermodynamic potential 
i7 is a simple BCS sum over different bands. In Section III 
we consider the problem of paramagnetic pair-breaking 
in multiband superconductors and show that a new gap- 
less state is energetically stable in low magnetic fields 
in some region of model parameters. We investigate the 
stability and the magnetic properties of this new gapless 
state in an s-wave two-band superconductor analytically 
at T — 0, and provide the details for the B ~ T phase 
diagram and the low temperature critical point that sep- 
arates the new partially gapless state from fully gapped 
state. In Section IV we present our conclusions. 



II. THE MULTIBAND MODEL IN THE 
WEAK-COUPLING LIMIT. 

In a standard BCS approach [73l| the pairing interac- 
tion can always be eliminated in favor of a single en- 
ergy scale, Tc, giving rise to the well-known universal- 
ity of the BCS theory. The critical temperature Tc is 
the only parameter that determines thermodynamic, ki- 
netic, and other properties in the weak coupling limit. 
For example, the superconducting gap at zero temper- 
ature, A(0), is related to Tc in s-wave superconductors 
by the universal law, A(0) = {tt/^)Tc ~ 1.76Tc. Similar 
universality is applicable to non s-wave superconductors 
as well, even though the universal ratio A{0) /Tc depends 
on the type of pairing. For d-wave pairing, the weak 
coupling ratio of the maximum gap amplitude at T = 
to Tc is A{0) /Tc = 2'KTc/jy/^ ~ 2.U^. Since a num- 



ber of interaction parameters are involved in a multi- 
band model, there is no such universal relation between 
Aq,(0) on different Fermi surfaces and T^- Thus, Tc can- 
not be the only parameter that describes the properties 
of multiband superconductors. Geilikman, Zaitsev, and 
KresinjSl, |3| showed using the method of Pokrovskii[75J 
that some universality is left in the weak coupling multi- 
band model. First, the physical properties, such as ther- 
modynamics, are often additive over different bands. Sec- 
ond, the ratios of the gap amplitudes on different Fermi 
surfaces are temperature-independent in the weak cou- 
pling limit. Third, Geilikman et aL[3| found that all 
physical properties of a BCS multiband superconductor 



can be expressed in terms of the transition temperature 
Tc, and other quantities measurable in the normal state, 
such as the ratios of densities of states on different Fermi 
surfaces, and the temperature-independent ratio gap am- 
plitudes. The gap amplitudes themselves, however, are 
not universal. In this section we introduce the multi- 
band model, review some of the results of Geilikman 
et al.[Z\ that we will use in other sections, and general- 
ize their weak-coupling solution to the case of arbitrary 
anisotropic pairing. 

The Hamiltonian for several separate Fermi surfaces 
has the following form (see, for example, Ref. [72J ) : 



Hci = ^e(fe)aj,^(fc)a„„(fe) + -^ ^ Vapa,^Ak,k')al^^{-k)ai^,^{k)apa-Ak')ap<TA-k') 



(1) 



fcjfc' a0(Ti 



r 



where ai-4 are spin indices, Vaisai^iil^,^') corresponds 
to the model interaction for a pair of electrons from band 
/3 with quasimomentum k' to band a with quasimomen- 
tum k. The latter can be written in the following form: 



Va/3CTi_4(fc,fe') = Val3{k, k'){icry)crj^a2{i(^v)l 



(2) 



Note that the unbalanced pairing terms [53 ISJ, llMl are 
not present in this model BCS formulation, since the 
mismatch between different Fermi surfaces is usually too 
large (of the order of ^ eV) for these terms to be relevant. 
We consider below any possible type of the superconduct- 
ing state. As usual, the interactions in the model BCS 
Hamiltonian are taken in a factorized form. 



Vo^pik-k') = x{^,e)v^pxy 



(3) 



where x((/9, 9) is the appropriate irreducible representa- 
tion, normalized to unity: 



in 



x{^,oW 



1. 



For example, for the s-wave pairing 

while for the d-wave pairing, 

xi'f^O) = \/2cos(2(/j). 



(4) 



(5) 



(6) 



The energy gaps on each Fermi surface Ac(fc) can 
be easily expressed in terms of the interaction ma- 
trix Vaf3{k;k') and the anomalous Gor'kov functions 

Ff3{iujn,k): 

A„(fc) = -T J2 V^p{k,k')Fp{iUr,,k') (7) 



Fi3{iuJn,k) 



Ap{k) 



■e + \MkW 



(8) 



The transition temperature Tc is determined by the lin- 
earized gap equation Eq.([7]) for the gap amplitude 



AA^,e)^Ac,xi^,e), 



Aa = ^ Xa/lAp In 



nTn 



where 



Aa/3 = —VaflVfS, 



(9) 



(10) 



(11) 



while the eigenvector A/3 corresponding to the largest 
eigenvalue of the interaction matrix Aq./3 determines the 
ratios between gaps on different Fermi surfaces set at Tc- 
Here i/p are the densities of states (DOS) per one spin 
direction for various Fermi surfaces. While the system 
of gap equations Eq.([7]) seems to give a temperature- 
dependent ratio for gaps on different Fermi surfaces, this 
result would violate the BCS logarithmic approximation. 
Within the logarithmic accuracy, the ratios between the 
gaps on different Fermi surfaces is set at Tc from Eq. pO)) , 
and is a temperature-independent constant set by the 
interaction matrix Eq. pTj) . the highest eigenvalue gap 
eigenvector [3|, \^ . Then the multiband problem can be 
parameterized in terms of Tc and the gap ratios, simi- 
lar to how this is done in the single-band BCS model. 
An additional related difhculty comes from the fact that 
the kernel for the system of Fredholm integral equations 
Eq.© is asymmetric, i.e., Xa/s ^ Xfia for Vot^ vp. Gener- 
alizing the approach of Pokrovskii[75| and following Gei- 
likman et a/.[3j], we introduce new variables and a sym- 
metric kernel yb^p = /i/3a, 



Ca 



VaVal3-/i^- 



(12) 
(13) 



As usual, the gap equation in terms of the universal scale 
Tc can be obtained by subtracting Eq. pU]) from Eq.©. 



Writing the result in new variables, we find: 



/3n 



\xiv, 



1 



I 



u^?.+e 



13 



T 



(14) 



Let us now multiply Eg. p^ by C,a, and sum over a, using The gap equation is then considerably simplified, and can 
the symmetry of ^ap and Eq. pU)) . be written in the universal form: 



/ , CafJ'aP — 2_^ M^aCa 



C/3 



In^ 



(15) 



T 






d^ujj 



\xiv,d)\' 



lx(¥', 



In 



T 



(16) 



Here 



^/3 






(17) 



are constant ratios, determined by the gap eigenvector at 
Tc and the corresponding densities of states. The coeffi- 
cients Uj3 are automatically normalized: 



E^ 



= 1. 



(18) 



As demonstrated in the above derivation, all physi- 
cal properties in the multiband BCS model can be ex- 
pressed in terms of the transition temperature Tc, the 
relevant DOS, and temperature-independent gap ampli- 
tudes. The gap amplitudes themselves, however, are not 
universal. The ratios of the gaps on different Fermi sur- 
faces are non-universal constants, determined by the rele- 
vant interaction matrix (the highest-eigenvalue eigenvec- 
tor for Aa/3 in Eq. pO|) l. Eq.((T6]) can be easily solved at 
T = 0. Introducing the usual normalization for the gap. 



^BCS - 

and defining ta as 



TlTr 



-{|x(¥',l 



•in|x(v,e)|) 



AqO = ta^BCS, 



(19) 



(20) 



we find: 



E^ 



, Intn 



(21) 



which sets a constraint on the amplitude of the gap eigen- 
vector at T = 0. Ea. p9| gives a standard expression for 
the zero-temperature amplitude of the energy gap for any 
single-band BCS superconductor. For example. 



ttTc 
Ascs = — - ^ 1.76Tc, s-wave, (22) 

7 

V2ttTc 



^BCS 



jVe 



1.51Tc, d — wave. (23) 



Note that due to the normalization condition Eq. Q , our 
definition for the d-wave gap amplitude is different from 
the standard definition by a factor of \/2. The temper- 
ature dependence of the gap amplitude A(r) is deter- 
mined by the universal gap equation Eq. (|16p . For the 
general case multiband case it differs from the single- 
band BCS temperature dependence: 



In- 



T 
TV. 



EE^ 

a n— 



lx(^,^)l^ 



V(" + 0.5)2 + (i„A(r)/27rT)2|x(<^, ( 



n -1-0.5 



(24) 



Since the ratios of the gaps on different Fermi surfaces are 
set at Tc, one can obtain thermodynamic potential inte- 
grating the gap equation Eq. pl]) over the single coupling 
constant, Tr'. 



Thermodynamics is then given by a sum of standard 
weak coupling expressions for separate bands expressed 
in terms of temperature-dependent energy gaps Aq (T) : 



^S — ^N = 



dT, 



rpi y.^'^'-^a I rpi 



T' ^ 



T 



T' 



(25) 



ns~vi 



N 



-^TY,vJ^ul + ^i\x{v,ew 



^ui + ^l\x{v,ew 



2\uj„ 



(26) 



Integrating the above expression over uj at zero temper- 
ature, one finds the familiarfS'l factorized result for the 



ground state energy: 



Es-En=- 



v„I^l 



L^ 



(|x(¥',0)P)o = -A|csE 



VaK 



(27) 



r 



where we have used the normalization condition for 
x((/?, 0) given by Eq.(|3]). The multiband gap equation 
does produce overall change for the relevant quantities. 



such as, for example, specific heat jump at transition 
temperature, T^, obtained, for example, in Ref.[3j: 



J 



AC 
C 



12 



Eo^aAi 



12 (E„^aA2)^ 



7C(3)Ea^a (Ea^^A^) 7C(3)Ea^a 



Ea ^oA4 



(28) 



r 



The above result of Geilikman et al. 3\ is applicable to 
any anisotropic or unconventional superconductors be- 
longing to a one-dimensional representation of the point 
group [Ty, [73], which includes d-wave multiband super- 
conductors. For superconductors belonging to a degener- 
ate representation of the point group, the corresponding 
generalized formula is given by 



AC 



AC 



C 



1 



(Ec 



J^^A? 



BCS 



Ea ^" Ea ^o,Ai 



(29) 



where [AC/C]bcs is the weak coupling value of the spe- 
cific heat jump at Tc for a given multiband representa- 
tion, given, for example, by Kuznetsova and the author 
inRef.[7|. 

As a simple example, let us consider the two-band 
case. The expression for the constant ratio of the energy 
gaps on different Fermi surfaces is then very easily ob- 
tained from Eq. pO)) in terms of the interaction constants 
Eg. pT)) (also see, for example, Ref.[23j): 



A2(T) _ A2(T,) ^ 2Ai2 

Ai(T) Ai{Tc) A22 - All + V(^ii - ^22)2 + 4A12A21 

I 



(30) 



The other relevant parameter of the BCS model is uf, defined by Eq. (fT7|) : 



2 ^1 1 2 

i^l + i^2S^ 



(31) 



Since the gaps on all Fermi surfaces have the same tem- 
perature dependence, let us introduce a normalized gap 
amplitude A as given by Eq. ((20|) . 



A„(T) = i,A(T), A{T = 0) = Abcs- 



(32) 



We then easily find the expression for zero-temperature 
gaps on the two Fermi surfaces: 



t-\ = s 



U^s'' 



(33) 



The temperature dependence of the gaps on the two 
Fermi surfaces must be the same. In addition to Tc, 
thermodynamics is completely determined by two other 
parameters, s and ui, which can be easily found from 
experiment. 

We conclude that the universality of the weak-coupling 
BCS-like model is applicable to the multiband case. How- 
ever, some non-universal constants, which depend on 
the interactions, do enter the problem as temperature- 
independent parameters. The total number of indepen- 
dent parameters for an m-band BCS model is then sig- 
nificantly reduced to 2m — 1 measurable constants: Tc, 
771—1 independent constant DOS ratios Va/vp^ and ni—1 
independent constant gap ratios A^/A^. A significant 
temperature dependence of the energy gap ratios on dif- 
ferent Fermi surfaces as, for example, observed in MgB2, 
then corresponds to the failure of the weak coupling log- 
arithmic approximation. 



III. PARTIAL GAPLESS STATE IN MAGNETIC 
FIELD IN S-WAVE SUPERCONDUCTORS. 

In this section we consider the unbalanced pairing 
problem [4 ll. |42| for a multiband s-wave superconductor 



J 



placed in paramagnetic magnetic field. The author and 
Gor'kovJTlj have recently shown that a new stable gapless 
state appears in small magnetic fields. The emergence of 
this state for s-wave superconductors is accompanied by 
zero-temperature first order phase transition. The s-wave 
case can be solved analytically at T = 0. In what follows, 
we present our results for the stability of the gapless state 
in s-wave multiband superconductors and the details of 
the B ~ T phase diagram. 



A. The nature of the gapless state 



We start by considering the general solution for the 
BCS gap equation for a multiband superconductor placed 
in magnetic field. We neglect the diamagnetic effects 
{Hc2)- The presence of paramagnetic magnetic field re- 
sults in an additional Pauli term in the Hamiltonian for 
each seperate band: 



Hr, 



J 



{k){I- 



,(fe), I = tiBB. (34) 



'2 ^9cca'„ 

kija 

Here ga is the g- factor for each band a, taken to be 
isotropic in this section. The multiband gap equations 
in magnetic field have the same form Eq.([7|, with the 
Gor'kov function F given by 



^a(w„,p) 



-tayAkip) 



iiLUn-Ia^y~Uky~\A^ip) 



(35) 



r 



The diagonalization procedure of Pokrovskii[75J and Gei- 
likman et al. 3] is applicable in the presence of arbitrary 
magnetic field, i.e., the ratios of the gaps on different 
Fermi surfaces do not change as functions of both tem- 



J 



perature and magnetic field, and are determined by the 
gap eigenvector at Tc- For example. Eg. ([50]) for the two- 
band problem in the presence of magnetic field [Tlj can 
be written as: 



A2(T,B) A2(Te, 



Ai(T,B) Ai(T< 



cO) 



2Ai 



^^22 — Ml + V ('^ll ^ -^22) + 4A12A21 
I 



(36) 



where Tco ^Tc{B = 0). 



Thus, as it is the case with the energy gaps and ther- 
modynamics, it is possible to obtain a complete solu- 



tion of the multiband problem in magnetic field in terms 
of Tc, gap ratios, and DOS ratios. The linearized gap 



J 



equations in magnetic field produce the instability curve 
Tc{B), which takes the following form: 



ln^=*fi 






u^Re 



<_ 



r 



UnTr 



(37) 



where Ua are band-dependent constants given by Eq. ()17|) . 
Note that in the free electron case, ga = 2, the multiband 
instability curve Eq. ([57)1 is reduced to the familiar single- 
band result, in agreement with earlier results [7l|: 



In — ^ = * - 



Re 



, 1 I 
* hi 

2 2TTTr 



(38) 



The re-entrant behavior of T^ with increased magnetic 
field normally indicates the possibility of first-order phase 
transitions on the B — T phase diagram. In a single- 



band problem the homogeneous gap equation gives rise 
to a second unstable solution 59] in magnetic fields close 
to the paramagnetic limit, also known as the Sarma[43| 
state. The instability is resolved in favor of an inhomo- 
geneous L0 FFj4ll \4^ state. As we have shown in the 
recent letter [71[, in the two-band case one has three dif- 
ferent solutions. In the general rn-band case, the number 
of solutions will be m+1. The T = solution of the gap 
equation will change form at / = Aq(-B). For the s-wave 
case all solutions can be written out analytically as: 



J 




V(g^/)2-4A^(/)2+g,/ 

2AaQ 



r 



n 

a,gcI<2Af,{I) 



A^/) 



i/30 



= 1. 



(39) 



Here a and (3 are band indices, u^ are given by Eq. (fT7|) . 
while Aqo correspond to the solution of the multiband 
gap equation at T = without the magnetic field, 
Eq. ipO)) . Stability of these solutions can be inferred from 



the ground-state energy, which can be easily obtained 
by the integration of the gap equation over the coupling 
constant T^ as in Eq. ((25| : 



Es — Eno — — y ^ 



yaAJIf 



Y, vpgpl^igpiy - 4A^(/)2 

/3,g0/>2A3(/) 



(40) 



Here -Eato is the normal-state energy in the absence of 
magnetic field. 

We can also write out the solution of the multiband 
gap equation at finite temperatures, since there is only 
a single temperature and field-dependent variable, the 
same for all Aa{T,B): 



SiT,B) 



A„(T,B) 



Aa{T = 0,B = 0) 



(41) 



The gap equation then takes a simple form: 



lniSiT,B)]^j:ulfJ^^,'4) 



(42) 



where 



fo{x,y) = 



pOO 



dt 



Vt^ 



X 



x^ 



' sinh iV^ , ^ 

cosh (y) + cosh (V? 



,2) 



- X^ 



(43) 



which is the same function that appears in the single- 
band model. The expression for thermodynamic poten- 
tial in magnetic field is an analytic continuation of the 



J 



corresponding expression in zero field to iuJ„ 
{gal/ 2), which is also factorizable 



lUJn 
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with 
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fi{x,v) = ^ / t^dt I cosh" 



\/W+x^ ~ y 



cosh 



VW^Tx^ + y 
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Here 



ilN{I) = ^m ^ 7 X! '^aidal)^ 



(46) 



is the normal state energy in exchange magnetic field /. 
From the general analysis of Eq. (|1(I)) one can see that 
the solution with gal > 2Aa{I) for all a is always unsta- 
ble, similar to the unstable Sarma state of the single-band 
problem^]. The low-field solution, gal < 2Aq,(/) = 
2Aao for all a, on the other hand, corresponds the multi- 
band BCS solution at / = 0, which is stable in low enough 
magnetic fields. The other solutions of the gap equation 
correspond to partial Sarma states, illustrated on Fig[T] 
for the two-band case and ga = 2. The energy spectrum 



of the system for excitations near each FS is given by the 
poles of the Green's function Gai^^n, l^)- 



1 



Ea (fc) = VCa(fe)2-f|A„(fc)|2 + - gala' 



(47) 



As the magnetic field exceeds the value of the small- 
est gap, a strip of unpaired fully polarized quasiparticles 
forms in the vicinity of the corresponding Fermi surface, 
giving rise to a paramagnetic magnetic moment. Sim- 
ilar to the single-band problem, the LOFF state com- 
petes with homogeneous solutions in high magnetic fields. 
A generalization of the instability curve Eq. (l37|) can be 
written for an arbitrary inhomogeneous q- vector as: 



Tco \2 



E^ 



Re 



* 



. ggl + 2{vFaq) 
' 47rT, 



(48) 



r 



where vpa is the Fermi velocity for band a. Tc{I) is then 
found as a maximum with respect to q. The high-field 
phase transition is always second order. As the field is 
lowered, the regions of stability of various phases and the 
exact nature of the LOFF state in the general case, when 
the energy gaps are of the same order of magnitude, can 
only be obtained numerically. 

The zero-temperature phase transition in partially 
gapless state in s-wave multiband superconductors is 
accompanied by the appearance of the paramagnetic 
moment |71j. In non s-wave multiband superconductors, 
however, magnetic moment is always present in low mag- 
netic fields due to the nodes in the energy spectrum. The 
transition to partial Sarma state then corresponds to a 
crossover from the nodal regime to a regime with a full 
open Fermi surface. In completely isotropic situation, the 
zero-temperature phase transition to partial Sarma state 



in s-wave superconductors is always 1-st order, corre- 
sponding to the appearance of a finite paramagnetic mag- 
netic moment. Effects of spin-orbit, non-spherical Fermi 
surface, or gap anisotropy can turn this phase transition 
into a smooth crossover. Similarly, due to effects of gap 
anisotropy, the first-order zero-temperature phase transi- 
tion in multiband d-wave superconductors exists only in 
a certain region of parameters jTlj. The first-order phase 
transition, if present, disappears above a certain critical 
temperature, Tcr- 

The analysis of the B ^ T phase diagram for multi- 
band superconductors depends on the same number of 
parameters as the analysis of other properties considered 
in previous section. Thus, for example, energetic stability 
of various gapless state will depend on the value of these 
additional parameters. We thus consider for simplicity a 
two-band model, which has only two additional physical 
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FIG. 1: Multiband partial Sarnia state at T = 0, charac- 
terized by fully polarized unpaired electrons near the Fermi 
surface of the driven band. The state is stabilized by the 
presence of the gap on the primary Fermi surface. 
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constants, the ratio of the densities of states on two Fermi 
surfaces, 1^2/ '^i, and the ratio of the energy gaps A2/A1. 
We also assume for simplicity gi = g2 — 2 and only 
consider homogeneous solutions, since the LOFF state 
will depend on the shape of both Fermi surfaces. Fig. 
[2^ shows an example B — T phase diagram for homo- 
geneous phases for the two-band isotropic s-wave model 
in exchange field and parameters (A2/A1) — 0.2 and 
(1^2/^1) = 25 obtained with the help of the above expres- 
sions for the temperature-dependent gaps and thermody- 
namic potential Jig. The instability line in Fig. [2^ for the 
transition into a uniform partially gapless superconduct- 
ing state has the characteristic reversal behavior, which 
indicates the presence of either a LOFF state [4l|,|42| or a 
first order phase transition. The LOFF instability line for 
two equal circular Fermi surfaces with different masses in 
two dimensions is shown by the red dashed line. When 
only the uniform superconducting states are considered, 
the the partially gapless state will be separated from the 
normal state in high magnetic field and low temperatures 
by a first order Clogston-typc transition. At low fields it 
is separated from the fully gapped superconducting state 
by a first order line that ends in a critical point. Fig. 
?? shows the details of the first order phase transition 
from fully gapped superconducting state to partially gap- 
less superconducting state. The first order transition is 



FIG. 2; (a) B — T phase diagram for a 2D two-band supercon- 
ductor with two equal circular Fermi surfaces (ppi = PF2) and 
(A2/A1) = 0.2, (m-2.lm\) = 25. The blue line corresponds to 
the first order phase transition into the partial Sarma state 
from the normal state (right side, the Clogston limit) or fully 
gapped superconducting state (left side). The LOFF insta- 
bility is shown by a dashed red line. The instability for the 
normal/uniform superconducting state is shown by the black 
line, (b) Details of low temperature first order transition from 
fully gapped to gapless superconducting state. The black 
line marks the first order phase transition. The red dashed 
lines correspond to the boundaries of the region of metastable 
states, where the hysteretic behavior is expected to occur. 



shown by the solid black line that ends in a critical point 
at T = TcR '^ O.OSSTco. The dashed red fines mark the 
boundaries metastable region in the vicinity of the first 
order phase transition where unstable solutions of the 
gap equation are present. 

Fig. [3] shows the solution of the multiband gap equa- 
tion Eq. (j42p as a function of temperature and magnetic 
field. The unusual reversal behavior of the energy gap re- 
flects the presence of a first order phase transition. While 
there are three different solutions at a given field in a re- 
gion near this phase transition, only one of these solutions 
is stable. Such a region corresponds to the region of the 
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FIG. 3: Numerical solution of the two-band BCS gap equation 
A(T, B) at different temperatures for parameters in Fig[2] 
Three solutions for the gap equation at low fields and tem- 
peratures reflect the existence of unstable region near the flrst 
order phase transitions from fully gapped state into partially 
gapless state. As the temperature is raised above the criti- 
cal temperature, the first order region disappears, and there 
is only one solution for the gap equation in low fields. The 
high-field behavior corresponds to the usual first order unsta- 
ble Clogston limit. 



FIG. 4: Free energy Qs of the superconducting state for a 
two-band superconductor in magnetic field at different tem- 
peratures for the parameters in Fig[2l The reversal behavior 
at low temperatures T < Tcr = O.OSSTco indicates the pres- 
ence of the first order phase transition from the fully gapped 
superconducting state into a partially gapless gapless super- 
conducting state. The first order transition from the par- 
tially gapped superconducting state to normal state in high 
magnetic fields happens at the point when the condensation 
energy f7s(r,/)=0. 



first order phase transition where the hysteresis exists. 
The stable solution corresponds to the minimum of the 
Free energy, Eq. (P^ . shown in Fig. [H The shape of the 
Gibbs Free energy at T < Tcr = O.OSSTco also reflects 
the presence of the metastable states near the first order 
phase transition into the gapless superconducting state. 
The system always picks the lowest energy. Thus, the free 
energy of the system as a function of magnetic field has a 
slope change at the point of first order phase transition, 
where the energy gap and paramagnetic magnetization 
have a corresponding jump. 

The zero-temperature phase transition into the new 
state is characterized by a metamagnetic jump of mag- 
netization. At T = the magnetic moment appears 
sharply, from A/ = to M ^ 0. At a finite temperature 
small magnetization is present due to thermal popula- 
tion of the band with smaller energy gap. The first order 
phase transition line disappears at T > Tcr, as the quasi- 
particle states above the gapless state become thermally 
populated. Thus, one has two metamagnetic transitions 
shown in Fig. [31 one corresponding to the transition into 
partially gapless Sarma state, the other is the first order 
phase transition from the partially gapped Sarma state 
into the normal state. 

The analysis of the above equations can be done nu- 
merically for any parameters Vijvx and A2/A1. For the 
s-wave case, however, analytical results for the energy 
and the magnetic moment can be obtained at T = 0. In 
a special case when the critical point is located at low 
enough temperatures, the critical region is fully deter- 



mined by the second Fermi surface and depends only on 
one energy scale. We consider these simplified results 
below. 

Finally, we note that the partial gapless state is not 
always present on the B — T phase diagram; for some 
parameters of the two-band model, as shown below, the 
first order transition happens directly from the uniform 
superconducting state into the normal state. 



B. Magnetic properties of the partial gapless state 

at r = o 



In this section we consider the low field transition to 
partial gapless state. As we have seen in the previous sec- 
tion, the gap equation for the two-band model has three 
different solutions and depends on two constant param- 
eters, the ratio of the energy gaps on the two Fermi sur- 
faces s, determined only by interactions, Eq. (|36p . and the 
ratio of the densities of states for the two bands. How- 
ever, analytical expressions are more conveniently writ- 
ten in terms of s and a combination of these parameters 
ul, defined in Eq. ([3T|) . 

The uniform low-field solution at T = is the same as 
the solution without the magnetic field. The first order 
phase transition appears near I ~ A20 . It is thus conve- 
nient to rewrite the solution in terms of A20 and other 
parameters of the second Fermi surface. For example, 
the energy of the superconducting state at T = and 
5 = is 
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FIG. 5: Magnetic moment normalized to magnetic moment in 
the normal state at different temperatures for the parameters 
in Fig. [21 The first order jump in magnetization from zero to 
finite magnetic moment in the lower fields and temperatures 
T < Tcr = 0.058Tco corresponds to a transition into partially 
gapless state. The sharp transition is absent at temperatures 
above the critical temperature, as quasiparticle states above 
the second smaller gap become thermally populated. 



where the solution changes abruptly from A = 1 to the 
smaller gap branch of Eq. (|53p . The magnetization in 
the partial Sarma state as a function of A is also easily 
determined, since it is also just a fraction of the normal 
state magnetization on the second Fermi surface: 



M 



-2/iBi/2^tanh(u2 In A) 



(56) 



Eqs. (j53p and ((56)) determine the field dependence of 
magnetization in the partial gapless state parametrically. 
Magnetization has a jump from to a finite value at the 
first order phase transition Acr- The transition point 
can be easily found analytically in a particular case of 
w| <C 1, corresponding to "induced" superconductivity 
on the second Fermi surfacefTlj. Indeed, introducing 



TA = A - 1, 



(57) 



Eg. ([53]) transforms to: 






(58) 





EsQ „ , ^2A2o. 


(49) 


Introducing 






A 


A2(S,T) Ai(B,r) 

A20 Aio 
/ 

A20 


(50) 
(51) 



we find a simple expression for / in the partially gapped 
state: 



J=-A(A"2 +A""2 



/ = A cosh (^2^ In A). 
The expression for the energy also simplifies: 



-^ = A + 2uliVP - A2 
Eso 



(52) 
(53) 

(54) 



The 1-st order transition point is determined from the 
smallest A < 1 solution of the transcendental equation, 

1 - ul sinh (2u^2 In Acr) = A^^^ (55) 



while for the condensation energy we easily find the fol- 
lowing expression: 



Es o 4 o 

Eso "" 3ii| ^' 



(59) 



The cubic terms in the energy lead to a first order phase 
transition 



III at 



Tic 



-3ut/8, 



(60) 



where the energy gap ta changes abruptly from ta = 
to TA = — 3m|/2. The magnetization changes abruptly 
at Tjcr from zero to 



Mcr = -U2^JLBV2A20■ 



(61) 



High-field Clogston limit and energetic stability 
of the gapless state at T = 0. 



We now turn to the first order transition in high mag- 
netic field, or the modification for the Clogston criterion 
for the two-band model. The standard Clogston Crite- 
rion involves a transition from the uniform superconduct- 
ing state to the normal state, determined by 
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EN{I)^-{yi+l^2)lllog 



(62) 
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FIG. 6; The region of parameters for the two-band model for 
which the T = gapless phase is present. 



We thus obtain 



^clog 



1 /z^iA|o + J^A|o 



1^1 + V2 



(63) 



The transition first order transition from gapped two- 
band superconductor to the normal state only happens 
at / = Iciog in the absence of the partially gapless state. 
Since the partially gapless state is different from the uni- 
form gapped state, the high field transition to the normal 
state will also happen in a different magnetic field, now 
determined by 



Es{I) = En{I) 
After simple calculations, we find: 



clogl 



UiS 



2 



Ul(l-s2) + Vui-2s2 ' 



(64) 



(65) 



where ui = -^/l — u^, and the new Clogston limit Idogi 
then determined by Eq. ([53]) that connects A and / in the 
partially gapless state. 

It is obvious that the condition 



^cloql ^ ^c 



(66) 



where Icr is the magnetic field for the T = transition 
from uniform into partially gapless state, must be met 
for the partially gapless state to be present on the phase 
diagram. In particular, it is obvious that for 



s > 



Ml 

V2 



(67) 



the partially gapless state is definitely not present. The 
accurate condition is given by Eq. ([55|) . which upon sub- 
stitution of Eq. ([S5)) into Eq. ([55)1 takes the following form: 



1 



; sinh 



(2^2^ In Ac 



■)>A- 



(68) 



The region of parameters of the two-band model where 
the partially gapless state exists given by Eq. ([55)) is 
shown in Fig. [Sj We see that in order for the gapless 
state to be present and easily observable, the second gap 
has to be much smaller than the first gap, while the den- 
sity of states on the second Fermi surface should be quite 
a bit larger than on the first one, otherwise the entropy 
change associated with the first order transition will be 
very small and not easily detectable. Such transition will 
also be turned into a crossover in the presence of even 
small gap anisotropy, impurities, non-zero temperature, 
or other effects |7g|. 

We have not considered the details of the high-field 
stripe LOFFJill, \^ region beyond the LOFF instabil- 
ity in the normal state given by Eq. (j48p and shown in 
Fig. [21 since this calculation depends on the details of 
the shape of the Fermi surfaces, their dimensionality, and 
thus a number of additional parameters [80|- Unlike the 
Clogston limit, the LOFF region is very non- universal 
and sensitive to defects. We note, however, that the 
LOFF stripe phase will also be unusual, since it forms in 
high enough magnetic fields in the vicinity of the first or- 
der transition between partially gapless state and the nor- 
mal state. In particular, the inhomogeneous LOFF state 
will likely involve superconducting stripe order made out 
of the partially gapless state. It is also not completely 
clear whether the phase transitions between the partially 
gapless superconducting state and the LOFF state and 
between the LOFF state and the normal state will be 
first or second order. 



D. The low- field critical region. 

The mathematics of the low-field critical region near 
the T = first order phase transition into the partially 
gapless superconducting state and the associated ther- 
modynamics is rather bulky, and can only be studied 
numerically for the general case (see Figsl5][5]). Never- 
theless, as it was shown in the previous subsections, only 
parameters for the second Fermi surface are relevant for 
the first order phase transition at T = 0. Unfortunately, 
both Fermi surfaces determine the critical region at fi- 
nite temperatures. However, when ^2 <C 1, the critical 
region lies at very low temperatures, and thus only the 
quasiparticle excitations near the gapless or nearly gap- 
less state on the second Fermi surface are important. The 
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results for the critical region then depend only on param- 
eters of the second Fermi surface only, and thus it can be 
found in a universal form. We note that the special case 
1*2 <?C 1, corresponding to superconductivity driven by a 
single Fermi surface and induced by interactions on other 
Fermi surfaces is, perhaps, most common for multiband 
superconductors. In this section we consider such a weak 
first order transition, assuming that both Fermi surfaces 
are completely isotropic. 

The T = first order phase transition in this limit 
happens near / ~ A20 • To describe this phase transition 
and thermodynamics near it, following Ref. [Tlj, it is 
convenient to introduce new dimensionless variables that 
correspond to deviation of the second energy gap and the 
magnetic field from A20: 



I- A 



Tl 



20 



TA 



_ A2(T,/)-A2 



i20 



and dimensionless temperature 

T 

" A^' 



(69) 



(70) 



(71) 



The solution for the metamagnetic transition at T = 
in this limit is given by Eqs. (|58l) . ([59|) . It is not difhcult 
to extend this solution to finite temperatures. The gap 
equation then takes the following form: 



TA = -W2 



'V~t 



^Jxdx 



cosh^ [: 



"^+"^2^ 



(72) 



Thermodynamics near the critical point is simply given 
by the integral of the gap equation, 

^s - ^so - 4|f]so| / ' T^{Ti)dTi. (73) 



The critical point can be found by differentiating Eg. ((72 
and solving the following two equations: 



drA dT\ 



0. 



The critical temperature is given by: 

with 

_ f"" ^tanh(x + A2)dx 

1 ^ / 2 ' 

Jo cosh {x + A-i) 

where A^ is the solution of 
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(74) 
(75) 
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^ l-2sinh"(a; + A2) 
cosh (X + A2) 



The energy gap and magnetic field at the critical point 
are given by 



'^A,, 



~Axu\ 



\fxdx 



cosh {x + A2) 



Solving for Ai and yl2 numerically, we obtain: 



(78) 
(79) 



tcr 


= 0.3129m|, 


(80) 


TA,cr 


= -0.7541ut 


(81) 


Tl,cr 


= -0.40714- 


(82) 



It is not difhcult to see that the solutions for the 
gap equation, the energy, and the magnetic moment in 
the critical region depend on single energy scale Tcr, or 
m|A2o. Thus, universal numerical results for the critical 
region can be obtained. 

In Fig. [7] we show the first order transition line on 
the B — T phase diagram in universal units u^ r/ and 
T/(A2ou|). Fig. [8]shows the solution of the gap equation 
in the critical region at different temperatures in univer- 
sal units U2 Tj and U2 ta. The gap equation has three 
different solutions. The black line is the line of first order 
phase transitions of Fig. [7] in these coordinates. The two 
Ta on this line at a given field tj correspond to the two 
coexisting superconducting states on the first order line 
at a given temperature. The thermodynamic potential 
ft = (fls/flso) — 1 is shown in reduced units Om^^ in 
Fig. [51 The behavior of the gap as a function of mag- 
netic field results, as usual, in the presence of metastable 
states in the vicinity of the first order phase transition. 
Finally, Fig llOl shows the jump in paramagnetic magne- 
tization at different temperatures below the critical tem- 
perature Tcr in reduced units of M/{^b'^2A.2oU2)- As 
the temperature is raised to Tcr, the first order jump 
disappears. 



IV. CONCLUSIONS. 

Let us now briefiy summarize our main results. We 
have performed a detailed calculation of energetically 
stable homogeneous superconducting states in the un- 
balanced pairing problem for s-wave multiband super- 
conductors. Our analysis shows that this problem differs 
from one for single-band superconductors, and that quali- 
tatively new partially gapless states may be present in the 
low temperature region of the B — T phase diagram. The 
new states are characterized by a gapless Fermi spectrum, 
open Fermi surfaces, and a finite paramagnetic magnetic 
moment. The phase transition between fully gapped and 
gapless superconducting states in magnetic field at T = 
is a metamagnetic first order phase transition, which cor- 
responds to a sharp jump in magnetization on one of the 
Fermi surfaces that becomes gapless. 

The superconducting order is present on both Fermi 
surfaces, as shown in Fig. [1] At finite low temperatures 
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FIG. 7: The T — B phase diagram near the first order tran- 
sition to partially gapless Sarma state in the special case of 
weak first order transition, marked by a red solid line. 



FIG. 9: Gibbs free energy of the superconducting state near 
the first order phase transition to the partially gapless state 
for different temperatures T < Tcr- The unusual behavior of 
the Gibbs potential indicates the presence of metastable states 
near the first order transition from fully gapped to partially 
gapless superconducting state. 
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FIG. 8: Solution for the gap equation in the low field critical 
region as a function of magnetic field and temperature. The 
black line is the first order transition line for different tem- 
peratures in B — A coordinates, also shown in Fig[7]in B — T 
coordinates. 




FIG. 10: First order metamagnetic transition into partially 
gapless state at different temperatures T < Tcr- The transi- 
tion disappears and becomes a crossover at T > Tcr- 



the metamagnetic first transition results a first order line 
on the B — T phase diagram that ends in a critical point, 
as shown in FiglH The presence of the gapless super- 
conductivity also modifies the high-field Clogston limit. 
While the new state is analogous to the one studied in 
Ref. [55], it is not the same. Unlike the situation encoun- 
tered in Bose condensation or high energy physics, unbal- 
anced pairing of different species of fermions in supercon- 
ductors is energetically unfavorable because of the large 
difference of the corresponding Fermi surfaces. Pairing in 
multiband superconductors is associated with each Fermi 
surface separately, although the gaps on different Fermi 
surfaces are related by the interaction. In the weak cou- 
pling logarithmic scheme that ratio is temperature- and 



field-independent. Nevertheless, we found that, similar 
to Ref . [55] , the gapless state is most visible when the su- 
perconducting gap on the heavier band is driven by the 
superconducting transition on the lighter band. 

Strong anisotropy of the B — T diagram, which in- 
dicates a quasi-2D nature of 2H-NbSe2, and its multi- 
gap supcrconductivityf23, [U, [23 makes a case for a 
possible first order phase transition to a ground state 
of this type in low magnetic fields parallel to the plane 
in this material 7l| . An unusual first order phase tran- 
sition has been indeed observed in this material [8l| for 
thermal conductivity measurements in low magnetic field 
H ^ IQkOe ^ Hc2 parallel to the basal plane. This 
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first order phase transition is inconsistent with explana- 
tions involving vortex lattice melting [Slj. The magnetic 
field at which this transition occurs is consistent with the 
value of the small energy gap ^ 0.1 — 0.2meV observed 
in the photoemission experiments [21| . The low- field first- 
order phase transition in NbSe2 was found to be strongly 
anisotropic, and the hysteretic behavior of thermal con- 
ductivity disappears for certain field directions 81], a be- 
havior not expected in a simple multiband model consid- 
ered above, where a first order phase transition occurs 
for all field directions in the basal plane. However, in 
the presence of CDW the CI symmetry is broken [82], 
which itself has been shown to lead t o a strong in-plane 
anisotropy of the B—T phase diagram[!36]. Similar effects 
must be present for the low-field metamagnetic phase 
transition as well. Gap anisotropy, g-factor anisotropy, 
and spin-orbit coupling tend to wipe out the first order 
line on the B — T phase diagram, and turn it in a smooth 
crossover |7l|. We have recently found that the order of 
the phase transition in the presence of spin orbit interac- 
tion may, too, be dependent on direction. As a result of 
two different terms in the spin orbit interaction, the first 
order phase transition becomes very anisotropic. It is 
present for some field directions and turns into a smooth 
crossover for other field directions [73]. 

Experimental observation of the new gapless state is 
subject to the usual difficulties associated with the ob- 
servation of paramagnetic pair breaking and the LOFF 



state in superconductors. Namely, the orbital effects 
leading to i?c2 are almost always present, even in strongly 
quasi-2D materials in magnetic fields parallel to the 2D 
planes. Perhaps, an ideal realization of this state would 
be surface superconductivity or superconductivity in thin 
films in fields parallel to the surface. We note, however, 
that unlike the LOFF state, which is difficult to observe, 
since it quickly disappears in the presence of impurities 
or orbital effects, the new homogeneous partially gapless 
state is more robust, and will appear in many multiband 
strongly quasi-2D s-wave superconductors in the mixed 
state as well, provided that the upper critical field is close 
enough to the Clogston limit. The details of such first 
order transition in the mixed state will be similar to the 
physics considered above. Measurements of the specific 
heat in applied field are the most direct way to observe 
the low-field first order phase transition in s-wave multi- 
band superconductors, such as 2II-NbSe9 [2]| . 
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